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Pulsars are good clocks in the universe. One fundamental question is that why they are good clocks? This is related to the braking 
mechanism of pulsars. Nowadays pulsar timing is done with unprecedented accuracy. More pulsars have braking indices measured. 
The period derivative of intermittent pulsars and magnetars can vary by a factor of several. However, during pulsar studies, the 
magnetic dipole braking in vacuum is still often assumed. It is shown that the fundamental assumption of magnetic dipole braking 
(vacuum condition) does not exist and it is not consistent with the observations. The physical torque must consider the presence of 
the pulsar magnetosphere. Among various efforts, the wind braking model can explain many observations of pulsars and magnetars 
in a unified way. It is also consistent with the up-to-date observations. It is time for a paradigm shift in pulsar studies: from magnetic 
dipole braking to wind braking. As one alternative to the magnetospheric model, the fallback disk model is also discussed. 
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1 Introduction 

Pulsars are rotating magnetized neutron stars. They are good 
clocks in the universe. Because of their high precision tim¬ 
ing, pulsars are good probes of interstellar medium and mag¬ 
netic field ll29l . and good probes of gravitational wave back¬ 
ground ED. At the same time, pulsars have many other ap¬ 
plications, e.g. pulsars as clocks etc. Among these applica¬ 
tions, one fundamental question is that why pulsars are good 
clocks? This is related to the emission mechanism of pul¬ 
sar multi-wave radiations and the spin down of pulsars. For 
the spin down mechanism of pulsars, one may ask, taking the 
Crab pulsar as an example, “how to explain its period deriva¬ 
tive, braking index etc?” The period and period derivative of 
various pulsar-like objects are shown in figure [D 

Pulsars must have magnetospheres (a system of particles 
surrounding the central neutron star). In the magnetosphere, 
there is particle acceleration. The radiation of these parti¬ 
cles is responsible for the pulsar emissions. However, the 
magnetic dipole braking in vacuum is often assumed in pul¬ 


sar timing studies, which is one dilemma in current pulsar 
studies. In the following, it is shown that the fundamental 
assumption of magnetic dipole braking (vacuum condition) 
does not exist, and it is not consistent with the pulsar ob¬ 
servations. A physical spin down mechanism must be based 
on the existence of a magnetosphere ll25l . One candidate at 
present is the wind braking model 08011931 . 

Section two is about the early understanding of pulsars. 
Observational progresses are provided in section three, in¬ 
cluding the braking index of pulsars, the spin down behaviors 
of intermittent pulsars and magnetars. Theoretical progresses 
are given in section four, including pseudo-magnetic dipole 
braking (various modifications of magnetic dipole braking), 
wind braking of pulsars, wind braking of magnetars and 
some discussions of the fallback disk model. Summary and 
prospects are presented in section five. The following ma¬ 
terials (especially section two) are influenced by three text¬ 
books Il541l721l94ll . During the preparation of this paper, we 
read the work of ||8](mainly chapter two there). One aim of 
this paper and |[8l is the same: trying to prove that the mag- 
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netic dipole braking in vacuum is incomplete (even irrele- are used in this paper. The corresponding rotational energy 
vant) for the study of pulsar spin down. loss rate is 



Figure 1 Period and period derivative diagram of pulsars, including nor¬ 
mal pulsars (black points), magnetars (blue squares, empty squares for radio 
loud magnetars), X-ray dim isolated neutron stars (green diamonds), central 
compact objects (light blue circles), rotating radio transients (red stars) and 
intermittent pulsars (magenta triangles). Updated from Figure 1 in [83] 


2 A bite of pulsar early researches 

When Hewish et al. Il^ first reported their discovery of the 
pulsating radio source (i.e. pulsar), the oscillations of white 
dwarfs or neutron stars were thought to be responsible for the 
observations. It was Gold Il2^ who first pointed out the ro¬ 
tating neutron stars as the origin of pulsars (the rotating neu¬ 
tron star scenario was not welcomed by the pulsar commu¬ 
nity when it was first proposed). If the radio emission comes 
from the rotational energy of the neutron star, then two natu¬ 
ral predictions will be made ||2^ : (1) pulsars will slow down 
gradually; (2) more rapidly rotating pulsars are expected (the 
pulsar period at that time ranges from 0.25 s to 1.33 s). Later 
discovery of the Vela pulsar (with period of 89 ms), the Crab 
pulsar (with period 33 ms), especially the slowdown rate of 
the Crab pulsar (about one part in 2400 per year) has con¬ 
firmed the rotating neutron star origin of pulsars Il24l . 

With the period and period derivative of the Crab pulsar, 
the central star’s rotational energy is 

£,ot = i/Q2 = 27t2/U (1) 

- 2xl04®^erg=«2xl04®erg, (2) 

where I is the neutron star’s moment of inertia (which is or¬ 
der of 10^^ g cm^, for a star mass 1-2 Mq and radius about 
10 km), Q. is the angular velocity, and P is the rotation period. 
Following the tradition of theoretical astrophysics, CGS units 


IFrotl = -£™t=47t2/^ (3) 

= 4x lO^^^^ergs^'a: 5 X lO^^ergs^*. (4) 

It is much higher than the radio luminosity. The energy out¬ 
put meets the energy budget required by the Crab nebula. 
Therefore, from the beginning of pulsar astronomy, it is al¬ 
ready known that the radiation energy is only a small frac¬ 
tion of the rotational energy loss rate (especially for the radio 
emissions). Most of the rotational energy is taken away by 
the particles and electromagnetic waves. This mixture of par¬ 
ticles and waves is named as particle wind (or wind) in the 
following. This particle wind may be visible in the form of 
pulsar wind nebulae Il22ll . 


2.1 Rotating magnetic dipole in vacuum 


No mathematical formula was presented in Gold’s theoreti¬ 
cal paper (this is Gold’s style). The magnetic dipole radiation 
formula was employed by Pacini to consider the extraction of 
rotational energy and slowdown of neutron stars before and 
after the discovery of pulsars. In 1967, just before the discov¬ 
ery of pulsars, Pacini proposed the idea of extracting ro¬ 
tational energy of neutron stars by magnetic dipole radiation. 
This may provide a continuous energy supply required by the 
Crab nebula observations. For a rotating magnetic dipole in 
vacuum, the radiation frequency is equal to the angular speed 
of the neutron star. The corresponding radiation luminosity 
is 

2u2q4 

Ea = ^^3 sin^ a, (5) 

where ^ is dipole magnetic moment, a is the angle between 
the rotational axis and the magnetic axis (i.e. the inclination 
angle), and c is the speed of light. The dipole magnetic mo¬ 
ment is related to the equatorial surface magnetic field and 
polar surface magnetic field Bp as: fi - BJi^ - 1 f2BpR^ {R 
is the neutron star radius). The magnetic field at the magnetic 
pole is two times as that at the magnetic equator Bp = 2Bs (for 
a dipole geometry). Equation Q can be rewritten in terms of 
magnetic field strength 


Ea 


2B}R^Q^ 


3c3 


sin^ a 


BlR^Pf 

6c4 


sin^ a. 


( 6 ) 


The dipole magnetic moment is about 10^° G cm^ for a mag¬ 
netic field of 10'^ G. The corresponding dipole radiation lu¬ 
minosity is about 10^®-10"^° erg s“' for an angular speed of 
10 ^-10^ rads-'. 

After the discovery of pulsars, ll^ considered the mag¬ 
netic dipole braking of pulsars further. By equating equation 
(O and (Is), the period evolution of the neutron star is 


P = 


87t^//^sin^a 
37^3 P 


(7) 
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If the timing parameters are known, the surface magnetic field 
can be obtained from equation (|2l) 

BsSina = 3.2xlO'®V^G. (8) 

Assuming an inclination angle of 90 degrees, the commonly 
employed expression for the characteristic magnetic field of 
pulsars is obtained 

Bc(equator) = 3.2 x lO'® '/pPG. (9) 

Therefore, the characteristic magnetic field is the equatorial 
magnetic field by equaling all the braking torque to an orthog¬ 
onal rotator in vacuum. Furthermore, for the study of pulsar 
emissions, the magnetic field strength at the magnetic pole 
is more relevant. It is two times larger than the equatorial 
surface magnetic field 

Bc(pole) = 6.4 X 10^'^ V^G. (10) 

The rotational evolution of a pulsar can be obtained by in¬ 
tegrating equation Q 


Pit) 



IbTT^/i^ sin^ a 

31? 


(t - to) 


+ 6.16x 10 tait - to), 


( 11 ) 

( 12 ) 


where Pq is the initial rotational period at time to, and Z?s,i 2 is 
the surface equatorial magnetic field in units of lO'^ G. Equa¬ 
tion (fT2l i is obtained for (f - to) in units of years. For a given 
set of parameters [Po, to, B^, a), when t is small, the rota¬ 
tional period P(t) changes very little P(t) Pq. During the 
late time, when t is large, the period increases with the time 
as P(t) oc fb2 -pjjg evolution of period derivative, rotational 
energy, and rotational energy loss rate etc can all be obtained 
directly from equation (fT2l i. From equation (fTTT i. the age of a 
pulsar is (using equation©) 

7’ = f-f0-Te(l-(y) j, (13) 


three effects will result in a braking index larger than three 
(see eauation(l29l)). However, the magnetic dipole braking 
assumes a rotating dipole in vacuum. This may make it not 
only incomplete but also irrelevant for the braking of real pul¬ 
sars ifTm . 


2.2 The presence of pulsar magnetospheres 


There are many pioneering works since ||25]| . The scale height 
of the neutron star atmosphere is about 1 cm. The gravita¬ 
tional binding energy of a proton (or electron) on the neu¬ 
tron star surface is about 100 MeV (or 0.1 MeV). However, 
the electromagnetic force will dominate over the gravita¬ 
tional force. Dimensional analysis shows that ESI, the ra¬ 
tio of electric force to the gravitational force for a proton 
is ~ 10®, where e is the absolute value of electron 

charge, G is the gravitational constant, M is the neutron star 
mass, and m is the proton mass. Therefore, a magnetosphere 
must be created around the central neutron star |I8]|25]. The 
magnetosphere tends to corotate with the central neutron star. 
However, the rotational velocity cannot exceed the speed of 
light. The light cylinder radius is defined as where the rota¬ 
tional velocity equals the speed of light 

7?,, = T = =4.8x 10®Pcm. (15) 

Q. 27T 

For the magnetic field lines closed within the light cylinder, 
they will corotate with the neutron star (“closed field line re¬ 
gions”). In the closed field line regions, in a steady state, 
the Forentz force vanishes EH (which follows the treatment 

of El) 

= (16) 

c 

where E is the electric field. The corresponding space charge 
density is (named as the Goldreich-Julian charge density) 




■E^ — 


n-B 


2nc 


The corresponding particle density is 

lie-1 X cm“^. 


(17) 

( 18 ) 


where 

P 

Tc = — (14) 

2P 

is defined as the characteristic age. If the magnetic dipole 
braking in vacuum is valid, and the initial period is much less 
than the present period Po P, then T - ll28l . On the 
other hand, if the initial period is close to the present period 
Po ~ P, then T «c Tc, i.e. the characteristic age may be much 
larger than the pulsar’s true age. The central compact objects 
may correspond to this case ll26ll . 

The magnetic dipole braking can be further developed by 
considering the presence of gravitational wave radiations and 
free precession due to deformation of the neutron star by 
its magnetic field, the decay of the magnetic field, and ra¬ 
diation reaction on the inclination angle IflTlIb^ . All these 


where is magnetic field component parallel to the rota¬ 
tional axis. The Goldreich-Julian density is the charge den¬ 
sity, which is the net value between positive and negative 
charges. The particle density of positive and negative charges 
can be much higher. 

However, this electric field (eg.lfTbll) is perpendicular to 
the magnetic field E • B - 0. In the corotating frame (an 
inertia frame corotating with the neutron star at some instan¬ 
taneous time), the electric field is - E + -xB - 0. In 
the closed field line regions, particles move the along field 
lines. The perpendicular electric field has no effect on parti¬ 
cle acceleration. The electric field responsible for the particle 
acceleration is the parallel component E • B 0 (which is a 
Forentz invariant). In the open field line regions, when the 
charge density deviates from the Goldreich-Julian density, it 
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will cause the appearance of a parallel electric field. Sepa¬ 
rating the total electric field in the open field line regions as 
E - E^cc + Eqj, where V • Eci = 47TpGj, the accelerating 
electric field satisfies the equation 

V-Eacc =47T(p-pGj)- (19) 


For a rotating dipole, the field line equation is r = rj sin^ 6, 
where rj is the maximum radial extension. When the max¬ 
imum radial extension equals the light cylinder radius, it is 
defined as the last closed field line (or the last open field line). 
The corresponding colatitude at the neutron star surface is the 
polar cap angle 


Maybe, the most simple case would be a gap, where there is 
no net charge, i.e. p -0. Near the neutron star surface, when 
the acceleration region is small, ea.(fT9b can be rewritten in a 
one-dimensional form 

d£|| 

^ =47t(p-pgj), (20) 

di 

where E\\ is the electric field component parallel to the mag¬ 
netic field line, s is the distance measured along the magnetic 
field line, and in the gap p = 0. The gap bottom (the neutron 
star surface) corresponds to s = 0. Denote the gap height as 
s - H. If the electric field at the top of the gap is defined 
as £’||(s - H) - 0, then the analytical solution for the elec¬ 
tric field can be obtained for the aligned rotator (the magnetic 
axis and the rotational axis are parallel to each other) ll9^ 
(adopting a similar procedure as GO)) 

2 QBp 

£|l =- -is-H). (21) 

c 

The potential difference between the top and bottom of the 
gap is 

A(f> = - -H^. (22) 

c 

For a pulsar with a rotational period of one second. Bp = 
10'^ G, and a gap height of H - 10"^ cm, the potential differ¬ 
ence is about Aip - 6 x lO'^ V. The gap height is determined 
by the complicated physical processes in the magnetosphere. 
It is related to the specific particle acceleration model. Var¬ 
ious physically motivated acceleration models are proposed 
since iTTOl . 

The magnetospheric torque can be estimated in several 
ways. The electromagnetic field far away from the neutron 
star may be of a transverse wave form 11^17211 . The local 
electric field will be the same as the magnetic field E ~ B (in 
CGS units). The corresponding Poynting flux is S ~ 

At the light cylinder, the energy loss rate is 

E ~ S47xRi = P , . (23) 

The dipole magnetic field decreases with the radius as B(r) 
Bp(Rlr)^. This relation is used in the deduction of the above 
equation. Equation (l23t has the same form as the vacuum 
dipole case (eq.®). 

Another way of calculating the magnetospheric torque is 
by counting the energy carried away by each outflowing par¬ 
ticle. It is determined by the geometry and particle acceler¬ 
ation in the open field line regions. The footprint of all the 
open field lines makes up a region, dubbed as the polar cap. 


6 »pc = sin' 




p-l/2^ 


(24) 


where the neutron star radius is taken as R = 10® cm. The 
polar cap radius is Rp^ = ROp^ ^ 1.4 x lO^^F *^^cm. For the 
electric field satisfying eg. lfTbl l. the corresponding potential 
drop between the polar cap edge and the magnetic pole is (it 
is actually the maximum acceleration potential for an aligned 
rotating dipole fTOl ) 


R^flBp ~ n 12 

AO^ax = - 6.6 X 10*2-^ V, (25) 

where Bp ,12 is the magnetic field in units of 10'^ G (in later 
sections, it will be simply written as B 12 ). The potential is 
expressed in units of volts. During numerical calculations, it 
should be converted back to CGS units. The particle flux in 
the polar cap region can be written as pc - kpq^c, where 
the charge density is parametrized as k times the Goldreich- 
Julian density. When k - \ and the acceleration potential 
equals the maximum acceleration potential, the correspond¬ 
ing energy loss rate due to the particle outflow is 


E = 27tBp^pGJcACl)niax = ' ^^6) 

It is also similar to the vacuum dipole case (eq.®). 

In the magnetosphere, particles must attain some critical 
acceleration potential (about lO'^-lO'^ V) in order to gener¬ 
ate secondary particles cniiBi. For a constant acceleration 
potential A(f) =constant, the corresponding energy loss rate is 
(after some deduction similar to ea.(l26ll) 


E = 


BpR^Q? 

— -A0. 

c 


(27) 


The slow down of a pulsar can be expressed in a power law 
form 1721 

Q oc -G", (28) 

where n is called the braking index, the minus sign corre¬ 
sponds to the slow down of the pulsar. Observationally, the 
braking index of a pulsar is defined as l5^ 


Gfi _ vv _ 2 BP 

Q2 ~ y2 - p2 ■ 


(29) 


In pulsar studies, the angular velocity Q = 27t/P = 271 v, 
pulse frequency v = 1/P, and the rotational period P are 
frequently used in different conditions. Pulsar researchers 
should have no difficulty in transforming one expression to 
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another. For the dipole braking case, the braking index is 
n - 3. The braking index is n = 1 for the magnetospheric 
torque in the case of constant acceleration potential (eq.dZTb). 
Therefore, the magnetic dipole braking and the magneto- 
spheric torque are different on the second order effect, i.e. the 
braking index. This defines the condition when the magnetic 
dipole braking approximation can be used. 

3 Observational progresses 

3.1 Braking index measurement of pulsars 

For a rotating dipole in vacuum, the theoretical braking in¬ 
dex is exactly three, if the magnetic field, moment of inertial, 
and inclination angle are constant. ED correctly measured 
the braking index of the Crab pulsar ; n - 2.515 + 0.005. It 
is different from the prediction of magnetic dipole braking. 
This strengthens the belief that the rotating dipole in vacuum 
is incomplete in modelling the spin down of pulsars. A sec¬ 
ond source had its braking reported in ll57l . For PSR B1509- 
58, its braking index was n - 2.83 + 0.03. It is similar to 
the Crab pulsar. In the presence of a particle wind, a smaller 
braking index « = 1 is expected EqI (and previous discus¬ 
sions, eg.dZTli). In the presence of a pulsar magnetosphere, 
there may naturally be some kind of particle outflow. These 
two early observations also indicate that it is hard to measure 
the braking indices of pulsars. It requires the determination of 
period second derivatives. The period second derivatives of 
most of the sources are dominated by “noise” processes ED. 

There are eight pulsars with braking index measured to 
the end of 2014 15^ . The braking index now ranges from 
n = 0.9 + 0.2 for PSR JI734-3333 HD to « = 2.839 + 0.003 
for PSR B1509-58 ll48l . In addition to the increase of the 
number of sources, the current braking index observations 
also show something new. 

1. The smallest braking index is for PSR J1734-3333: 
n - 0.9 + 0.2 lfT9l . It is consistent with a wind domi¬ 
nated magnetospheric torque ll^ . 

2. There are hints for the evolution of pulsar braking in¬ 
dex II 20 I . The braking index may evolve from about 
three for young sources to about one for old sources. 
The evolution of pulsar braking index should be con¬ 
sidered in theoretical models. 

3. The braking index is smaller during recovery from 
glitches. For PSR J1846-0258, its braking index de¬ 
creased from n - 2.65 + 0.01 to n - 2.16 + 0.13 after 
glitches Eq). The braking index of the Crab pulsar is 
about n - 2.3 during a glitch active epsiode ll56l . This 
may be viewed as evidence for glitch induced magne¬ 
tospheric activities in normal pulsars 1901 . 

More observations are required to confirm these interesting 
aspects. 


3.2 Spin down of intermittent pulsars 

In the early age of pulsar studies, some pulsars were found to 
stop working for several periods, i.e. nulling 0. Intermittent 
pulsars may be viewed as a special kind of nulling pulsars. 
The hrst intermittent pulsar PSR B1931-1-24 was discovered 
long ago. l3^ reported that PSR B1931-1-24 showed two dis¬ 
tinct states. In the “on” state, the pulsar can be detected and 
viewed as an ordinary pulsar. In the “off” state, the pulsar 
is not detected (various upper limits for the radio flux are re¬ 
ported). More intriguingly, the pulsar spin down rate is larger 
in the on state than that in the off state. The ratio of spin 
down rate between the on and off state is Pon/PoS = 1.5 
This means that there may be some additional particle out¬ 
flows during the on state. This particle outflow is responsible 
for both the radiation of radio emissions and larger spin down 
rate (because these outflowing particles will take away addi¬ 
tional angular momentum of the central neutron star iHTl l. 
The intermittent pulsars can have accurate determination of 
spin down rate for both the on and off state because their time 
scale of on and off state is very large. For PSR B1931-I-24, its 
on state lasts about 5 days. The off state lasts about 30 days. 
The duty cycle of the on state is about 20% (the fractional 
time when the pulsar is in the on state). 

Up to the end of 2014, three intermittent pulsars are re¬ 
ported lfT3llMll5n. The latter two intermittent pulsars can 
be in the off state for about one year. The spin down ra¬ 
tio ranges from PonIPoff = 1-5 (for PSR B1931H-24 ED) 
to Pon/PofF = 2.5 (for PSR J1841-0500 O)- Timing and 
radiation are two major aspects of pulsars. Intermittent pul¬ 
sars clearly show the correlation between timing and radiative 
behaviors. This kind of correlation may exist in all pulsars 
(only different by the magnitude of variations) 15^ . At the 
same time, efforts are made toward finding similar spin down 
behaviours in nulling pulsars (e.g. a larger spin down rate 
when the pulsar is on ||9^ ). Furthermore, rotating radio tran¬ 
sients ll58l and magnetar pulsed radio emission cni were also 
discovered. Along with intermittent pulsars, they all belong 
to the increasing diversity of transient radio pulsars. More 
transient sources may await to be discovered in future sur¬ 
veys (good news for future survey project, e.g. ||33). Future 
more accurate timing and multi-wave observations may un¬ 
veil the nature of these transient radio pulsars (also helpful 
for “persistent” radio pulsars). 

3.3 Timing behaviors of magnetars 

More than twenty magnetar candidates are discovered ll^ 
(http://www.physics.mcgill.ca/~pulsar/magnetar/main.html). 
They form a separate class of pulsars (neutron stars pow¬ 
ered by their magnetic energy IfTSl I. Magnetars manifest 
themselves as anomalous X-ray pulsars and soft gamma-ray 
repeaters II36II84I . Though their number is small, each source 
has its own peculiarities. 

Radio loud magnetars These sources belong to the in- 
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creasing number of transient magnetars. The first tran¬ 
sient magnetar XTE J1810-197 was reported in 2004 
El. Later it was reported to show pulsed radio emis¬ 
sions no). This discovery of radio emitting magne¬ 
tar bridged the gap between magnetars and normal ra¬ 
dio pulsars (the majority of rotation-powered pulsars). 
It also provides new clues to the study of pulsar ra¬ 
dio emissions (spectra, transient, and why radio emis¬ 
sions at all). Later multiwave observations found a de¬ 
creasing spin down rate (by a factor of three; here the 
spin down rate refers to the period derivative) while the 
star’s radio and X-ray luminosity keep decreasing IfTTl . 
Similar decreasing spin down rate was also found in 
another radio emitting magnetar PSR J1622-4950 ll4l 
when the star’s radio and X-ray flux decrease with time 
El . A decreasing spin down rate during outburst is nat¬ 
urally expected in the wind braking model of magne¬ 
tars ROl . 

For the radio magnetar IE 1547.0-5408, its spin down 
rate increases when its X-ray flux is decreasing with 
time 02. Similar negative correlation was also found 
in the radio magnetar near the Galactic center SGR 
J1745-2900. Its period derivative may have increased 
by a factor of two while the star’s X-ray luminosity is 
decreasing with time ESl . This negative correlation 
between radiative and timing behavior is contrary to 
the most simple idea of a magnetospheric torque. It 
may due to geometrical reasons ll85l or accumulation 
of magnetic energy in the magnetosphere (which may 
leads to outbursts in the future E3)- 

Low magnetic field magnetars Traditionally, magnetars 
were thought to be young neutron stars with a strong 
magnetic field Il40l . The difference between the dipole 
field and the multipole field is not clearly distinguished. 
This situation changed with the discovery of a magne¬ 
tar with dipole magnetic field lower than 7.5 x 10'^ G, 
i.e. a low magnetic field magnetar ED. The presence 
of a low magnetic field magnetar may be due to phys¬ 
ical IISTlI or geometrical reasons 117911 . It clearly shows 
that a strong dipole field is not the essential ingredient 
of magnetars. In the magnetar model, there are var¬ 
ious kinds of magnetic fields (core poloidal/toroidal 
field, crustal poloidal/toroidal field, magnetospheric 
poloidal/toroidal field, and local magnetic field do¬ 
mains etc). There are many parameters in the mag¬ 
netar model instead of one. More low magnetic field 
magnetars were discovered later II68II98I . 

Enhanced spin down rate and anti-glitch Many magne¬ 
tars show enhanced spin down rate during burst ac¬ 
tive episode (IE 2259-H586 ED; SGR 1806-20 J92ll l. 
Marginal evidence of net spin down of the central neu¬ 
tron star was found in SGR 1900-1-14 previously ED. 
Later clear evidence of net spin down was found in the 


continued monitoring of IE 2259-t-586 (dubbed as anti¬ 
glitch ID). Similar net spin down was also observed 
in PSR J1846-0258 (which showed some magnetar 
activities at that time Il49l ). Quasi-periodic torque vari¬ 
ation was found in IE 1048.1-5937 El- Anti-glitch 
may just be a period of enhanced spin down rate If82ll . 
The enhanced spin down rate may result from more ac¬ 
tivities in the magnetosphere. However, it is not clear 
at present why the magnetospheric activities can repeat 
quasi-periodically. 

4 Theoretical progresses 

4.1 Pseudo-magnetic dipole hraking 

In the magnetic dipole braking model, the rotational energy 
loss rate is proportional to sin^ a. The inclination angle will 
evolve to a smaller value El- The braking index will be 
larger than three considering the evolution of magnetic incli¬ 
nation angle (inconsistent with pulsar observations). In order 
to explain the observations (e.g. the braking index), various 
modifications of the magnetic dipole braking model are con¬ 
sidered. 

Generally, the spin down power law can be written as m 

Cl = -/f(f)Q^ (30) 

where the coefficient K{t) may dependent on the magnetic 
field, inclination angle, angular velocity etc. The correspond¬ 
ing braking index will be 

« = 3-—, (31) 

Tk 

where Tc is the characteristic age, and = Kjlk is defined as 
the characteristic variation time scale of K{t). If k > 0, then 

> 0. This will result in n < 3. This means that an increas¬ 
ing K will result in a braking index smaller than three. For the 
vacuum dipole case, k<x. ^ sin^ all. Therefore, an increasing 
magnetic dipole field strength lfT9ll . or an increasing inclina¬ 
tion angle ESll . or a decreasing moment of inertia 1971 will 
result in a braking index smaller than three. Furthermore, the 
increase of magnetic field and increase of inclination angle 
are coupled El ■ 

There are also other calculations done in the presence of 
a magnetosphere. At the same time, they can be under¬ 
stood (at least qualitatively) in the vacuum dipole case. In 
the presence of a pulsar magnetosphere, the rotational energy 
loss rate is made up of two components (at least mathemati¬ 
cally msiiiiiiii). They are the magnetic dipole component 
(which is proportional to sin^ a) and the particle wind com¬ 
ponent (different models have different expressions for this 
component). As the pulsar approaches the death line, the par¬ 
ticle wind component will gradually cease. For pulsars near 
the death line, their rotational energy loss rates are mainly 
dominated by the magnetic field perpendicular to the rota¬ 
tional axis, i.e. Bp sin a. For a small inclination angle, the real 
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surface magnetic field Bp can be much larger than the char¬ 
acteristic magnetic field, which is about Bp sin a. For the low 
magnetic field magnetar SGR 0418H-5729, its characteristic 
magnetic field is only about 7.5 x 10'^ G Wi\ . It already lies 
below the death line. For an inclination angle of 5 degrees, 
its surface dipole field can be as high 10*"^ G fTOl . Therefore, 
it can be a normal magnetar instead of a low magnetic field 
magnetar if its inclination angle is small. 

Magnetohydrodynamical simulations of pulsar magneto¬ 
sphere found a rotational energy loss rate proportional to 
(1 H- sin^ a) ITT^ . The term which depends on the inclination 
angle is the same as the magnetic dipole braking case, pro¬ 
portional to sin^ a. The evolution of inclination angle tends 
to reduce the rotational energy loss rate Il 66 ll . Then the incli¬ 
nation angle will also evolve to smaller value in the magne¬ 
tohydrodynamical simulation. The consequent braking index 
will also be larger than three. 

The above two cases are just two examples. Along with 
various modifications of magnetic dipole braking, they may 
be dubbed as “pseudo-magnetic dipole braking”. In the lit¬ 
erature, there are many theoretical and observational papers 
where the calculation or discussion are based on the magnetic 
dipole braking assumption. It may give the reader the illusion 
that the magnetic dipole braking model is well established or 
at least not too bad. In reality, it is not the case. The magnetic 
dipole braking model may be viewed as a pedagogical model, 
especially for newcomers in pulsar astronomy. 


4.2 Wind braking of pulsars 

Pulsars are oblique rotators in general. For a finite inclination 
angle, the magnetic moment can be decomposed into the per¬ 
pendicular component (perpendicular to the rotational axis) 
and the parallel component. The rotational energy loss rate 
due to the perpendicular magnetic moment may be approx¬ 
imated by the magnetic dipole radiation (see eq.Q, which 
may actually be some kind of particle outflow) 


£d 




sin^ a. 


The parallel component may be responsible for the particle 
acceleration no). The corresponding rotational energy loss 
due to the particle outflow is 19^ 

, Ad, 

Ep = InrpCpAcf, = ^^3—, (32) 


where Acp is the particle acceleration potential (depending on 
the specific acceleration model), and AO is the maximum ac¬ 
celeration potential. The total rotational energy loss rate may 
be the combination of the perpendicular component and the 
parallel component ||9^ : 


2//2q4 


E = 


2 2 I 

sin O'H-3-cos a\ — 

AO 


, 2 r >4 


2p^Q' 


(33) 


where rj = sin^ a + 3^ cos^ a. The magnetic dipole brak¬ 
ing in vacuum corresponds to rj - sin^ 0. The angular factor 
cos^ a appears by considering that the parallel component is 
mainly responsible for the particle acceleration. It is possi¬ 
ble that the factor cos^ a does not appear 1441 . However, the 
difference is only quantitative. Equation (l33l l is obtained by 
assuming that the particle density equals the Goldreich-Julian 
charge density and the effect of pulsar death is not considered. 
For the polar gap model CQl, the corresponding expression 
for T]is: T] - sin^ a + 4.96 x cos^ a Il44l . The 

corresponding expression can be obtained for each acceler¬ 
ation model. Therefore, precise pulsar timing observations 
may be used to distinguish between different particle accel¬ 
eration models. There are also other versions of magneto- 
spheric torque (similar to equation (l33T l. 181431 and references 
therein). Every model has its own merits. They can be used 
to check each other. 

Intermittent pulsars During the on state, the presence of 
some additional particles may be responsible for both 
the turn on of radio emission and the enhanced spin 
down rate. The magnetic dipole braking may be em¬ 
ployed to model the spin down during the off state. 
Then the ratio of spin down rate between the on and 
off state is mi; 

Qon 0 sin^ O'4-3(A(*/A<I)) cos^ a 
r = = —7— = -T- > 1. 

Goff sin a sin a 

(34) 

Therefore, in the pulsar wind model, the spin down 
rate of intermittent pulsars during the on state is always 
larger than that during the off state. Eor the vacuum gap 
acceleration model Col, the expression of 77 is known. 
Eor typical parameters, the corresponding ratio of spin 
down rate is a function of inclination angle, see Eig- 
ure|2] The predicted range of inclination angle can be 
checked by future observations. At the same time, the 
braking index of intermittent pulsar can also be pre¬ 
dicted. 

The Crab pulsar and other sources Equation ( 1^ may be 
further improved by considering the effect of particle 
density and pulsar death. Such a more perfect model 
of wind braking of pulsars is available in ll^ . It can 
be applied to the the Crab pulsar which has the most 
detailed timing observations. Eigure [3 shows the ro¬ 
tational evolution of the Crab pulsar in the wind brak¬ 
ing model. Many things can be seen in Eigure [3 Il39ll : 
(1) The spin down of the Crab pulsar will evolve from 
magnetic dipole braking dominated case (with brak¬ 
ing index close to three) to wind dominated case (with 
braking index about one). The observed value of brak¬ 
ing index (2.51 for the Crab pulsar) is the combined ef¬ 
fect of magnetic dipole braking and particle wind. Eor 
the whole pulsar population, their braking index can 
be naturally divided into two groups; one group has a 
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braking index close to three; the other group has a brak¬ 
ing index close to one. Pulsars in the second group will 
be older than those in the first group. (2) The Crab pul¬ 
sar (and all other pulsars) will not evolve to the cluster 
of magnetars. When the pulsar is sufficiently old, it will 
enter into the death valley. This is also true for the mag- 
netar population. Specific calculation for the low mag¬ 
netic field magnetar has already been done ll79l (using 
another version of wind braking). (3) When the parti¬ 
cle wind is stronger than the normal value (e.g. dur¬ 
ing glitches), the braking index will be smaller. This 
stronger particle will also contribute to some net spin 
down of the pulsar. This may explain the lower brak¬ 
ing index of the Crab pulsar during glitches. It may 
be viewed as glitches induced magnetospheric activi¬ 
ties in normal pulsars . (4) The wind braking model 
can also be applied to other sources, e.g. pulsars with 
braking index measured, and the magnetar population. 



Figure 2 Ratio of spin down rate for intermittent pulsars. The solid line is 
model calculation, for typical pai'ameters of: 5 = 5 x 10^^ G, P = Is. The 
dashed line is the observational range of spin down ratio, from 1.5 to 2.5. 
Adapted from figure 1 in E). 


■/□ magnetar • Rpp 
♦ XDINS ★RRAT 


▲ intermittent''pulsar ^ 



10 “’ 10 “’ 10 “' 10 “ lo' 

Period (s) 


Figure 3 Rotational evolution of the Crab pulsar in the wind braking model 
(red solid line). The eight sources with braking index measured are shown as 
grey squares with arrows marking their evolution direction. See figure^for 
explanations of various pulsar-like objects. Updated from figure 6 in (39) 

4.3 Wind braking of magnetars 

During the timing study of magnetars, the magnetic dipole 
braking in vacuum is also often employed. Furthermore, a 
pulsar is often claimed to be a magnetar when its character¬ 
istic magnetic field is higher than the quantum critical value 
(which is 4.4 x 10'^ G, when the electron cyclotron energy 
equals its rest mass energy). These claims do not consider the 
difference between magnetars and high magnetic field pulsars 
ED. Even if the putative star is indeed a magnetar, its charac¬ 
teristic magnetic field cannot be used as a proof of its magne¬ 
tar nature. This point is clearly demonstrated by the discovery 
of the low magnetic field magnetar SGR 0418-1-5729 (with 
a characteristic magnetic field less than 7.5 x lO'^G ill). 
Furthermore, there are other observations arguing against a 
strong dipole field in magnetars IHOll : the supernova ener¬ 
gies associated with magnetars are of normal value Il88l ; the 
non-detection of magnetars by Femu'-LAT II76II77I ; and the 
variable period derivative of magnetars (by a factor of several 
or more; on days or months time scale) etc. The variable spin 
down rate may be viewed as the direct evidence of magneto- 
spheric torques M75I80I . It is hard to image that the large scale 
dipole field can vary, especially on short time scale. While, 
the magnetospheric torque (e.g. a system of particle outflow) 
may change dramatically even on short time scale. The sys¬ 
tem of particles can change dramatically because they orig¬ 
inated from the magnetic energy. The X-ray luminosities of 
magnetars can vary significantly because it is also of mag¬ 
netic origin. Therefore, in the wind braking model of mag¬ 
netars, the timing and radiative behaviors of magnetars are 
correlated. This is commonly observed in magnetars 151 (and 
references therein). 
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The magnetic energy release may first be converted into a 
system of non-thermal particles. These particles are responsi¬ 
ble for the non-thermal radiation of magnetars. If the particle 
wind luminosity is Lp, then the rotational energy loss rate due 
to this particle wind is 

= (35, 

where Ed is the rotational energy loss rate due to a rotating 
dipole in vacuum. Two extremes of equation dTSl l are: 

Normal pulsar case The particle wind originates from the 
rotational energy Lp = -Emt- Then the rotational en¬ 
ergy loss rate due to the particle wind is order of mag¬ 
nitude the same as the vacuum dipole case. The effect 
of particle wind will mainly be seen in high order ef¬ 
fects, e.g. braking index ||9^ . Only in special sources, 
e.g. intermittent pulsars, the first order effect of particle 
wind can be seen El. 

Magnetar case The particle wind is from the the magnetic 
energy release, therefore it can be much higher than 
the rotational energy loss rate Lp » -Lrot- The rota¬ 
tional energy loss rate of magnetars will be dominated 
by the particle wind. A varying particle wind luminos¬ 
ity will result in a varying spin down rate of the puta¬ 
tive magnetar. A strong dipole field is not necessary in 
order to explain the timing observations of magnetars. 
In the wind braking scenario, magnetars are neutron 
stars with a strong multipole field. The strong multi¬ 
pole field is responsible for the persistent X-ray lumi¬ 
nosity, bursts (including giant flares), super-Eddington 
luminosity during bursts, and variable spin down etc. 

Therefore, equation ( l35l l represents a unification of wind 
braking of pulsars and wind braking of magnetars. Many con¬ 
sequences result from the wind braking of magnetars lISOl : 

1. The corresponding dipole magnetic field in the case of 
wind braking can be much lower (e.g. ten times lower 
or more) than the characteristic magnetic field. The 
surface dipole field at the magnetic pole is: 

Bp = 4.0x10^^- L^'{-^ G = 4.0x10'^^^^^— 

P p p,35 p/JOs P’35 

(36) 

where Lp 35 is the particle luminosity in units of 
1035 erg s"'. 

2. For magnetars with a low X-ray luminosity, they tend 
to have similar magnetospheres to those of normal pul¬ 
sars. Therefore, low luminosity magnetars are more 
likely to have radio emissions. 

3. There are two predictions of the wind braking model 
of magnetars. A magnetism-powered pulsar wind neb¬ 
ula is expected to surround the central magnetar. The 


braking index of magnetars will be smaller than three. 
These two predictions are both hard to varify observa- 
tionally. However, it is not impossible (especially in 
the future). 

Later observations are consistent with the wind braking 

model of magnetars. 

Swift J1822.3-1606 This is the second low magnetic field 
magnetar 1^ . However, different authors found dif¬ 
ferent period derivatives. While the observers were 
discussing whether this is caused by timing noise, 
IMl pointed out that the period derivative of Swift 
J1822.3-1606 may be decreasing with time. Accord¬ 
ing to equation (l35l l. for a short timescale of years, the 
star’s period derivative is related with the particle wind 
luminosity as: P oc . Therefore, a decreasing par¬ 
ticle luminosity after outburst will naturally result in a 
decreasing period derivative. IISTI predicted a long term 
period derivative of 1.9 x 10 (last paragraph, section 
2 there). Recent timing observations found a period 
derivative of (2.1 + 0.2) X 10“'"^ fTTlI . This is consistent 
with the wind braking prediction. 

SGR J1745-2900 This source lies near the Galactic centre 
and is radio loud Il69ll . During subsequent X-ray obser¬ 
vations, SGR J1745-2900 show an increasing period 
derivative (two times larger) when its X-ray luminosity 
keeps decreasing m- This negative correlation be¬ 
tween timing and radiation is hard to understand. Con¬ 
sidering the polar cap geometry of the particle wind, 
the particle wind may have a finite polar cap opening 
angle 0s- Then the corresponding period derivative will 
be proportional to P oc 0 /^ 3 ^ Therefore, a smaller po¬ 
lar cap angle will result in a larger period derivative. 
The change of polar cap opening angle will cause a 
negative correlation between X-ray luminosity and the 
spin down rate lf85l . In the wind braking model, SGR 
J1745-2900 will have a maximum spin down rate. Up 
to date timing observations are consistent with this pre¬ 
diction II16II86I . 

Anti-glitch The magnetar IE 2259-1-586 suffered a net 
spin down during an observational interval about two 
weeks. This net spin down is dubbed as anti-glitch lH. 
Previously, such spin down event also occurred in SGR 
1900-H14 m, and PSR J1846-0258 El. An en¬ 
hanced particle wind during the observational interval 
will take away some amount of additional angular mo¬ 
mentum. Therefore, the so call “anti-glitch” may just 
be a period of enhanced spin down of the central neu¬ 
tron star If82l . The X-ray flux is also higher than the 
quiescent level. This is consistent with the wind brak¬ 
ing interpretation of anti-glitch. 

Furthermore, the enhanced particle wind will result in 
a decreasing/varying braking index (if the star’s brak- 
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ing index can be observed). This kind of decreas¬ 
ing/varying braking index is indeed observed in the 
Crab pulsar Il39ll56ll90l . The opposite case of enhanced 
particle wind would be a period with a weaker particle 
outflow. This may correspond to the spin down behav¬ 
ior of intermittent pulsars Il44ll . Therefore, the wind 
braking model provides a unified explanation of (1) ro¬ 
tational evolution of the Crab pulsar, (2) spin down be¬ 
havior of intermittent pulsars, and (3) anti-glitch. 


4.4 Fallback disk model as an alternative to the magne- 
tospheric model 

The above discussions are mainly in the wind braking model 
(or magnetospheric model). For these peculiar timing behav¬ 
iors of pulsars and magnetars, there is always one alterna¬ 
tive; the fallback disk model. The fallback disk model is 
also commonly employed to explain various pulsar-like ob¬ 
jects II78II95I (as an alternative to the magnetar model). After 
the supernova explosion, some of the ejected material may try 
to fallback onto the neutron star. In the presence of some an¬ 
gular momentum, a disk around the central neutron star may 
be formed, i.e. fallback disk ll89l (and references therein). 
The observed pulsar braking index can be explained consider¬ 
ing the fallback disk torque 02111471 . The intermittent pulsar 
and other transient radio pulsars may result from the accre¬ 
tion of matter onto the neutron star from a fallback disk HD. 
However, quantitative calculations in the fallback disk model 
are still lacking (e.g. the spin down ratio of intermittent pul¬ 
sars has not been calculated quantitatively). More theoretical 
works are needed in order to confront with future multiwave 
observations. 

Anomalous X-ray pulsars and soft gamma-ray repeaters 
are proposed to be neutron star/fallback disk systems in the 
early age of magnetar research iniEi. Later timing observa¬ 
tions of magnetars can also be understood in the fallback disk 
model (“low magnetic field magnetar” ||2l; anti-glitch lITTII 
etc). The fallback disk model may account for the persis¬ 
tent X-ray emission and timing of anomalous X-ray pulsars 
and soft gamma-ray repeaters naturally (possibly more nat¬ 
ural than the magnetar model). In order to explain the gi¬ 
ant flares and the 10"^ times super-Eddington luminosity in 
the pulsating tail of giant flares, strong multipole fields as 
that of the magnetar model are required 1591 (and references 
therein). Considering that the central neutron star may be a 
quark star l95l . the quark star/fallback disk model may ex¬ 
plain both the persistent and burst properties of anomalous 
X-ray pulsars and soft gamma-ray repeaters EH). Future ob¬ 
servations (e.g. X-ray polarimetry 1521) may more clearly 
show which is closer to the truth, the magnetar model or the 
fallback disk model. 


5 Summary and prospects 

After about 50 years of pulsar discovery, pulsars timing are 
done with unprecedented accuracy. More pulsars have break¬ 
ing index measured, variation of period derivatives detected, 
timing noise characterised etc. The existence of intermittent 
pulsars directly shows the evidence of particle wind. Further¬ 
more, the timing behaviors of magnetars are always variable. 
However, the magnetic dipole braking in vacuum is still of¬ 
ten employed to explain the timing observations (not only by 
newcomers but also by experts). It is also well known that 
the fundamental assumption of dipole braking (vacuum con¬ 
dition) does not exist. Therefore, it is time for a paradigm 
shift in pulsar studies: from magnetic dipole braking to wind 
braking. The wind braking model is just one of the models 
for calculating the magnetospheric torque 0801I93II . The con¬ 
frontation between the wind braking model and pulsar timing 
observations is started; 

1. The braking of index of pulsars (generally 1 < n < 3) 
is the combined effect of magnetic dipole braking and 
particle wind. 

2. In the case of intermittent pulsars, the ratio of spin 
down rate between the on and off state is determined 
by the inclination angle. 

3. The effect of pulsar death is considered when mod¬ 
elling the long term rotational evolution of pulsars. 

More things are to be done in the future concerning the wind 
braking of pulsars: 

• The wind braking model can be applied to more 
sources, e.g. the other seven pulsars with braking in¬ 
dex measured, and the magnetar population. 

• The evolution of inclination angle in the presence of a 
particle wind. The evolution of rotational period and 
the inclination angle are coupled. The result may be 
compared with the Crab pulsar observations ll55ll . 

• The dependence of pulsar death line on the equation 
of state. The neutron star mass, radius and moment 
of inertia will affect the consequent dipole magnetic 
held, acceleration potential etc. The position of death 
line on the period and period-derivative diagram, and 
the rotational evolution will both depend on the pulsar 
equation of state. 

• The fluctuation of the magnetosphere may account for 
the timing noise of pulsars (e.g. HD). More detailed 
calculations should be done concerning the random 
fluctuation of the magnetosphere. 

• In the presence of a strong magnetic held, the neutron 
star will be somewhat deformed. This will introduce 
other torques, e.g. gravitational wave emissions. At 
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the same time, the deformed neutron star will also pre- 
cess (free or forced). During the life time of a pulsar, 
its magnetic field may also evolve, e.g. decay in the 
long run. Therefore, a thorough modelling should in¬ 
clude all these factors (some factor may only important 
in special cases). 

At present, the wind braking model of magnetars has 
achieved the aspects as follows; 

1. Unification of wind braking of pulsars and wind brak¬ 
ing of magnetars. 

2. A changing luminosity may explain the timing behav¬ 
ior of Swift J1822.3-1606. A changing geometry may 
correspond to the negative correlation between timing 
and X-ray luminosity in SGR J1745-2900. These are 
just two extreme cases. In reality, many factors can 
come into play. This may correspond to the various 
spin down behaviors of magnetars. 

3. The anti-glitch may be a period of enhanced spin down. 

The wind braking model of magnetars may be further devel¬ 
oped in the following aspects: 

• The particle luminosity is the main parameter in the 
wind braking model. At present, it is either assumed 
as lO^^ergs^' (or other constant values) or Lp = 
(which varies with time as that of the X-ray luminos¬ 
ity). In the future, the particle luminosity may be cal¬ 
culated for a given state of magnetosphere, e.g. the 
present attempt of 17] . Its evolution with time can also 
be calculated. This will enable the wind braking model 
to be applied to more sources. 

• The particle wind contributes both the radiation and 
braking torque of the magnetar. The present wind brak¬ 
ing model mainly concentrate on the braking torque. 
The radiation of these particles should also be ex¬ 
plored, e.g. radio, soft/hard X-ray. 

• The pulsar wind nebulae is due to the particle outflow 
in normal pulsars. In the case of magnetars, there may 
also be a wind nebula. This wind nebula may even 
powered by the magnetic energy of the magnetar. Both 
observation and theory are in need in this direction. 

The wind braking model is among the efforts to unify dif¬ 
ferent pulsar-like objects. The fallback disk model may pro¬ 
vide an alternative to the magnetospheric model. Accreting 
neutron stars form a distinct class of pulsars. Works for one 
special source are like blind men grabbing a piece of the ele¬ 
phant, i.e. the pulsar. The ultimate understanding of pulsars 
must involve a global view of all kinds of pulsars. 

The author would like to thank R. X. Xu for encouragement; F. F. Kou, L. 
Lin, and R. Yuen for reading the manuscript; and the Referee for comments 


about the fallback disk model. H.Tong is supported by the Xinjiang Bairen 
project, West Light Foundation of CAS (LHXZ201201), Qing Cu Hui of CAS, 
and 973 Program (2015CB857100). 

1 Alpar, M. A., 2001, ApJ, 554, 1245 

2 Alpar, M. A., Brian, U., Caliskan, S., 2011, ApJ, 732, L4 

3 Anderson, G. E., Gaensler, B. M., Slane, P. O., el al., 2012, ApJ, 751, 
53 

4 Archibald, R. R, Kaspi, V. M., Ng, C. Y., et al., 2013, Nature, 497, 591 

5 Archibald, R. R, Kaspi, V. M., Ng, C. Y, et al., 2015, ApJ, 800, 33 

6 Backer, D. C., 1970, Nature, 228, 42 

7 Beloborodov, A. M., 2009, ApJ, 703, 1044 

8 Beskin, V. S., 2010, MHD flows in compact astrophysical objects 
(Berlin: Springer) 

9 Bhattacharya, D., van den Heuvel, E. P. J., 1991, Physics Reports, 203, 
1 

10 Camilo, R, Ransom, S. M., Halpem, J. P, et al., 2006, Nature, 442, 892 

11 Camilo, R, Cognard, I., Ranson, S. M., et al., 2007, ApJ, 663, 497 

12 Camilo, R, Reynolds, J., Johnston, S., et al., 2008, ApJ, 679, 681 

13 Camilo, R, Ransom, S. M., Chatterjee, S., et al., 2012, ApJ, 746, 63 

14 Chatterjee, P, Hemquist, L., Narayan, R., 2000, ApJ, 534, 373 

15 Contopoulos, I., Spitkovsky, A., 2006, ApJ, 643, 1139 

16 Coti Zelti, R, Rea, N., Papitto, A., et al., 2015, MNRAS, 449, 2685 

17 Davis, L., Goldstein, M., 1970, ApJ, 159, L81 

18 Duncan, R. C., Thompson, C., 1992, ApJ, 392, L9 

19 Espinoza, C. M., Lyne, A. G., Kramer, M., et al., 2011, ApJ, 741, L13 

20 Espinoza, C. M., 2013, lAUS, 291, 195 

21 Eu, L., Li, X. D., 2013, ApJ, 775, 124 

22 Gaensler, B. M., Slane, P. O., 2006, ARA&A, 44, 17 

23 Gold, T., 1968, Nature, 218, 731 

24 Gold, T., 1969, Nature, 221, 25 

25 Goldreich, P, Julian, W. H., 1969, ApJ, 157, 869 

26 Gotthelf, E. V., Halpem, J. P., Alford, J., 2013, ApJ, 765, 58 

27 Groth, E. J., 1975, ApJS, 29, 453 

28 Gunn, J. E., Ostriker, J. P, 1969, Nature, 221, 454 

29 Han,J., 2013, lAUS, 291,223 

30 Hewish, A., Bell, S. J., Pilkington, J. D., et al., 1968, Nature, 217, 709 

31 Hobbs, G., Lyne, A. G., Kramer, M., 2010, MNRAS, 402, 1027 

32 Ibrahim, A. I., Mardwardt, C. B., Swank, J. H., et al., 2004, ApJ, 609, 
L21 

33 Karako-Argaman, C., Kaspi, V. M., Lynch, R. S., et al., 2015, ApJ, 
809, 67 

34 Kaspi, V. M., Gavriil, R P, Woods, P. M., et al., 2003, ApJ, 588, L93 

35 Kaspi, V. M., Archibald, R. R, Bhalerao, V., et al., 2014, ApJ, 786, 84 

36 Katz, J. L, 1982, ApJ, 260, 371 

37 Katz, J. L, 2014, Ap&SS, 349, 611 

38 Kramer, M., Lyne, A. G., O’Brien, J. T., et al., 2006, Science, 312, 549 

39 Kou, R R, Tong, H., 2015, MNRAS, 450, 1990 

40 Kouveliotou, C., Dieters, S., Strohmayer, T, et al., 1998, Nature, 393, 
235 

41 Lee, K.J., 2013, lAUS, 291, 189 

42 Levin, L., Bailes, M., Bates, S. D., et al., 2012, MNRAS, 422, 2489 

43 Li, J., Spitkovsky, A., Tchekhovsky, A., 2012, ApJ, 746, 60 

44 Li, L., Tong, H., Yan, W. M., et al., 2014, ApJ, 788, 16 

45 Li, X. D., 2006, ApJ, 646, L139 

46 Liu, X. W., Na, X. S., Xu, R. X., et al., 2011, ChPhL, 28, 019701 

47 Liu, X. W., Xu, R. X., Qiao, G. J., et al., 2014, RAA, 14, 85 

48 Livingstone, M. A., Kaspi, V. M., Gavriil, R P, et al., 2005, ApJ, 619, 
1046 

49 Livingstone, M. A., Kaspi, V. M., Gavriil, R P, 2010, ApJ, 710, 1710 


000000-12 


Tong H Sci China-Phys Mech Astron 


January (2015) Vol. 58 No. 1 


50 Livingstone, M. A., Ng, C. Y., Kaspi, V. M., et al., 2011, ApJ, 730, 66 

51 Lorimer, D. R., Lyne, A. G., McLaughlin, M. A., et al., 2012, ApJ, 758, 
141 

52 Lu, J. G., Xu, R. X., Feng, H., 2013, ChPhL, 30, 059501 

53 Lyne, A., Hobbs, G., Kramer, M., et al., 2010, Science, 329, 408 

54 Lyne, A., Graham-Smith, F., 2012, Pulsar astronomy (New York: Cam¬ 
bridge University Press) 

55 Lyne, A., Graham-Smith, F, Weltevrede, P, et al., 2013, Science, 342, 
598 

56 Lyne, A., Jordan, C. A., Graham-Smith, F., et al., 2015, MNRAS, 446, 
857 

57 Manchester, R. N., Durdin, J. M., Newton, L. M., 1985, Nature, 313, 
374 

58 McLaughlin, M. A., Lyne, A. G., Lorimer, D. R., et al., 2006, Nature, 
439, 817 

59 Mereghetti, S., 2008, Astronomy and Astrophysics Review, 15, 225 

60 Michel, F. C., 1969, ApJ, 158, 727 

61 Ng, C. Y, Kaspi, V. M., 2011, AIPC, 1379, 60 

62 Olausen, S. A., Kaspi, V. M., 2014, ApJS, 212, 6 

63 Ostriker, J. P, Gunn, J. E., 1969, ApJ, 157, 1395 

64 Pacini, F., 1967, Nature, 216, 567 

65 Pacini, F., 1968, Nature, 219, 145 

66 Philippov, A., Tchekhovsky, A., Li, J. G., 2014, MNRAS, 441, 1879 

67 Rea, N., Esposito, P, Turolla, R., et al., 2010, Science, 330, 944 

68 Rea, N., Israel, G. L., Esposito, P, et al., 2012, ApJ, 754, 27 

69 Rea, N., Esposito, P, Pons, J. A., et al., 2013, ApJ, 775, L34 

70 Ruderman, M. A., Sutherland, P. G., 1975, ApJ, 196, 51 

71 Scholz, P., Ka.spi, V. M., Gumming, A., 2014, ApJ, 786, 62 

72 Shapiro, S. L., Teukolsky, S. A., 1983, Black holes, white dwarfs, and 
neutron stars: the physics of compact objects (New York: Wiley) 

73 Spitkovsky, A., 2006, ApJ, 648, L51 


74 Sturrock, P. A., 1971, ApJ, 164, 529 

75 Thompson, C., Lyutikov, M., Kulkami, S. R., 2002, ApJ, 574, 332 

76 Tong, H., Song, L. M., Xu, R. X., 2010, ApJ, 725, L196 

77 Tong, H., Song, L. M., Xu, R. X., 2011, ApJ, 738, 31 

78 Tong, H., Xu, R. X., 2011, IJMPE, 20, 15 

79 Tong, H., Xu, R. X., 2012, ApJ, 757, LIO 

80 Tong, H., Xu, R. X., Song, L. M., Qiao, G. J., 2013, ApJ, 768, 144 

81 Tong, H., Xu, R. X., 2013, RAA, 13, 1207 

82 Tong, H., 2014, ApJ, 786, 86 

83 Tong, H., Wang, W., 2014, arXiv: 1406.6458 

84 Tong, H., Xu, R. X., 2014, AN, 335, 757 

85 Tong, H., 2015, RAA, 15, 1467 

86 Tong, H., 2015, arXiv: 1503.08904 

87 Turolla, R., Zane, S., Pons, J. A., et al., 2011, ApJ, 740, 105 

88 Vink, J., Kuiper, L., 2006, MNRAS, 370, L14 

89 Wang, Z. X., Chakrabarty, D., Kaplan, D. L., 2006, Nature, 440, 772 

90 Wang, J., Wang, N., Tong, H., et al., 2012, Ap&SS, 340, 307 

91 Woods, P. M., Koeveliotou, C., van Paradijs, J., et al., 1999, ApJ, 524, 
L55 

92 Woods, P. M., Kouveliotou, C., Finger, M. H., et al., 2007, ApJ, 654, 
470 

93 Xu, R. X., Qiao, G. J., 2001, ApJ, 561, L85 

94 Xu, R. X., 2006, An introduction to astrophysics (Beijing: Peking Uni¬ 
versity Press) 

95 Xu, R. X., 2007, Advances in Space Research, 40, 1453 

96 Young, N. J., Weltevrede. P, Stappers, B. W., et al., 2015, MNRAS, 
449, 1495 

97 Yue, Y. L., Xu, R. X., Zhu, W. W., 2007, Advances in Space Research, 
40, 1491 

98 Zhou, P, Chen, Y, Li, X. D., et al., 2014, ApJ, 781, L16 


